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The apparent low power in the CMB temperature anisotropy power spectrum derived from WMAP 
motivated us to consider the possibiUty of a non-trivial topology. We focus on simple spherical multi- 
connected manifolds and discuss their implications for the CMB in terms of the power spectrum, 
maps and the correlation matrix. We perform Bayesian model comparison against the fiducial best- 
fit ACDM based both on the power spectrum and the correlation matrix to assess their statistical 
significance. We find that the first year power spectrum shows a slight preference for the Truncated 
Cube space, but the 3-year data show no evidence for any of these spaces. 



INTRODUCTION 

The advent of high-precision cosmological observations 
has revived interest in the topology of the Universe. The 
latest measurements leading to questioning our assump- 
tion of a simply-connected topology came from WMAP 
[? ? ]. The anomalously low power at large scales in the 
temperature power spectrum, Ci, and suggestions of pre- 
ferred directions [? ? ] in the first year data could not be 
accounted for in the context of the standard flat ACDM 
model. A non-trivial topology could be a possible ex- 
planation. The observed low power, also reflected in the 
temperature correlation function, would arise naturally 
in this context as a consequence of the finite volume of 
the Universe. Thus, a number of (mainly flat), non-trivial 
topologies have been invoked, and various statistical tests 
employed to infer their validity (e.g. [??????? 
]). We study some of the simplest spherical non-trivial 
topologies. A closed geometry is marginally consistent 
with the data {Qk = I - = -0.02 ± 0.02). Moreover, 
the size of spherical manifolds (i.e., the length in various 
directions) is fixed with respect to the curvature radius, 
leaving only the curvature as a free parameter. Spher- 
ical manifolds have been discussed before, mainly with 
regards to their power spectra [? ? ? ]. Here, we expand 
upon this work, presenting full correlation matrices and 
maps and use Bayesian model comparison (first applied 
in cosmological model selection in [? ]) to assess their 
viability. 

We select our models using the geometrical degeneracy 
[? ], which states that universes of different spatial cur- 
vature have identical Q at small scales if their primordial 
power spectra, matter density and acoustic peak location 
parameter are the same. The degeneracy breaks at the 
large scales due to the ISW effect. We used this degen- 
eracy in [? ] to construct Q of simply-connected closed 
models showing reduced power at low £. We thus ensure 
our models exhibit the same high-^ Q as the WMAP best 
fit and focus on the low-^ regime, where topology plays 
an important role (we examine the 2 < £ < 10 region). 
The degeneracy constrains all cosmological parameters 
but the curvature, which is therefore related to the Hub- 
ble parameter; the two can be used as free parameters 



interchangeably [? ]. We explore -0.063 < < -0.008 
corresponding to 52 < i^o < 68 km/sec/Mpc. 

SPHERICAL NON-TRIVIAL TOPOLOGIES AND 
THEIR IMPRINTS 

We focus on the Quaternionic, Octahedral, Truncated 
Cube and Poincare spaces. Their topological proper- 
ties are listed in [? ? ]. Their fundamental domains 
are a four-sided prism, a regular octahedron, a trun- 
cated cube and a regular dodecahedron and the num- 
ber of fundamental domains tesselating the 3-sphere is 
8, 24, 48 and 120 respectively. These manifolds are 
rigid, indicating they have no other degrees of freedom 
aside from their orientation with respect to the coordi- 
nate system. The CMB temperature anisotropy spherical 
harmonic expansion coefficients and correlation matrix 

Cf^' {aimdl'rn') ^^^^ 

a,^ = ^A,(^)ypMEa,^e>, 
CL-' = ^A,(^)A,,(/3)P(/?)^^,^^^*,,^, 

where are the CMB transfer functions extracted from 
CMBFAST, P{l3) is the primordial power spectrum of 
perturbations and the e are Gaussian random numbers. 
The ^^^^ are particular to the space in question and en- 
code its topological properties. They are computed via 
a variant of the 'ghost method' [? ], which we will de- 
scribe in detail in [? ]. The /3 = (/c + 1)^/K are discrete 
wavenumbers where /c > 1 is an integer and K is the (un- 
normalized) curvature. Finally, we note that, while in a 
simply-connected trivial topology all > 2 correspond 
to eigenvalues of the Laplacian, this is not the case with 
the multi-connected topologies [? ]. In our analysis, we 
used < 41 for 2 < ^ < 10. 

Spherical multi-connected spaces exhibit reduced mode 
density (which can result in suppression of power) and 
anisotropic correlations at large scales. Fig. [T] shows the 
Ci of the four spaces of interest for Vtk = —0.017. 

A decrease in the size of the fundamental domain and 
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Figure 1: Power spectra of simply-/multi-connected spherical 
spaces for = —0.017 and a constant spectral index and 
from the first-year WMAP data. The drop at higher i is 
partly artificial and due to the small number of wavenumbers 
we used {/3 < 41). 



the number of allowed eigenmodes increases the sup- 
pression at large scales and the power spectrum appears 
more jagged. These effects become stronger if the curva- 
ture increases and at very high the Octahedral space 
shows a power deficit as well, supporting the argument 
[? ] that small well-proportioned universes suppress the 
quadrupole. The fit of the Q of the Truncated Cube 
and Poincare spaces is worth noting; however, it requires 
fine-tuning Qk- 

Our manifolds are intrinsically anisotropic, inducing 
apparent non-Gaussianity on CMB maps (i.e., if analyzed 
assuming isotropy the distribution of multipole moments 
at a single i might seem inconsistent with a univariate 
Gaussian). We constructed maps for 2 < ^ < 10 and 
show a sample in Fig. [2l There appears to be some 
structure in the maps of the two smallest spaces when 
= —0.063, but when the curvature decreases all maps 
appear interchangeable. Of course, we would have to 
analyze a large sample to identify any patterns. Upon 
decomposition of our realisations into individual multi- 
poles, we found occurrences of alignment; however, we 
need more realisations to draw any definite conclusions 
and the alignments we found may be isolated, patholog- 
ical features. 
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Figure 2: Maps for spherical non-trivial/trivial spaces. Top 
row: Qk = -0.063. Bottom row: Qk = -0.017. 



Signatures of multi-connectedness are much more ap- 
parent in the correlation matrix. Fig. [3] shows visuali- 
sations of a row of in pixel space. The Quater- 



nionic space shows no structure, the Octahedral space 
only shows correlations when Qq is high and the other two 
spaces always show correlations. The extent of these pat- 
terns grows as the fundamental domain becomes smaller 
and the topology more complex. These correlations ap- 
pear when the distance to the LSS, Xlss, becomes larger 
than the injectivity radius Rj (half the smallest geodesic 
between an object and its image) of a given topology. In 
our models, Xlss varies (in units of the curvature radius) 
from 0.649 (for Qk = -0.063) to 0.280 (for Qk = -0.007). 
The Ri for the Quaternionic, Octahedral, Truncated 
Cube and Poincare spaces are 0.785, 0.524, 0.397 and 
0.314 respectively [? ]. In our case, for the Quater- 
nionic space Ri > Xlss always, and for the Octahe- 
dral, Truncated Cube and Poincare spaces Rj > Xlss at 
^ -0.029, -0.012, -0.009 respectively. Beyond these 
limits, it is impossible to identify any of these topologies 
using the circles-in-the-sky method; however, they can 
be detected by our analysis, since off-diagonal terms still 
exist in the correlation matrix. 
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Figure 3: A row of Cf^ for spherical non-trivial/trivial 
spaces showing correlations with the point in the upper left. 
Top row: Qk = -0.063. Bottom row: Qk = -0.017. 

We also experimented with varying the spectral index 
n of the primordial power spectrum of perturbations: for 
a range of n from 0.5 to 1.2 and various values of Hq 
we found that a change in n tilts the power spectrum 
as expected, but leaves the patterns in the correlation 
matrix largely unchanged. 



MODEL COMPARISON 

The suppression of power at low £ suggests that these 
multi-connected manifolds could explain the WMAP 
findings. In order to assess their statistical significance, 
we perform Bayesian model comparison as in [? ]. 

Power Spectrum. We first compare the non-trivial 
topologies to the fiducial best fit model with respect to 
the 2 < £ < 10 part of the power spectrum. We use both 
the first- and three-year power spectra, because they dif- 
fer primarily as a result of a different analysis technique 
used by the WMAP team[? ]. Thus, we will be able to 
identify the impact of using different data analysis tech- 
niques on the evidence for spherical topologies. 

We express the likelihood as a function of the curvature 
and show it for both the first- and three-year WMAP 
data in Fig. [H In each case we have used the appropriate 
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likelihood code offered by the WMAP team (note that 
for the three-year data we use the low-^ component of 
the WMAP likelihood code). 
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Figure 4: Likelihood based on the power spectrum. 

We see that the first year data strongly favours a higher 
curvature (corresponding to lower values of Hq) and a 
smaller fundamental domain, as a smaller universe would 
show more suppression of power at low ^, in agreement 
with the data. In the case of the Poincare space, how- 
ever, the highly-closed regime ftk ^ —0.06 regime is not 
favoured, because the Q do not fit the data well. The 
same trend is observed in the likelihood based on the 
three-year data. We note that the drop in the likelihood 
at higher curvature is partly artificial and due to the fact 
that we only used /c < 40 in our calculations. 

We present the evidence for each topology in Table 
[H We impose a fiat prior on in the range [—0.063, 
-0.008] (52 km/sec/Mpc < Hq < 68 km/sec/Mpc). The 
upper bound is a result of computational limitations pre- 
venting us from getting accurate power spectra and cor- 
relation matrices at higher Qq. The lower bound is based 
on numerous suggestions that the Hubble constant is not 
likely to be any lower (see [? ] and references therein). 
The Octahedral and the Truncated Cube spaces are pre- 
ferred by the first year data with Bayes factors of 2 and 
8 respectively. The latter figure translates into definite, 
though not strong, evidence in favour of the Truncated 
Cube. However, the three-year data exclude all non- 
trivial topologies. 

If we switch from the agnostic fiat prior on 1]/^ to a 
Gaussian with Hq = 72 ± 10 km/sec/Mpc (allowing for 
systematic error) , we find a drop in the evidence for non- 
trivial topologies based on the first year data. This is 
expected, since this prior emphasizes higher values of i^o, 
where the likelihood is lower. Now only the Truncated 
Cube space is preferred, by a Bayes factor of just 2. When 
taking into account the three-year data a change in the 
prior does not have a significant effect. 

Correlation Matrix. Evidently, is still a free pa- 
rameter of our likelihood function; now, however, we 
must take into account the fact that Cf^ is not rotation- 
ally invariant. We can express the dependence on orien- 
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Table I: Top panel: Model comparison using Ci for the first 
(1st col.) and three-year (2nd col.) WMAP data. Bottom 
panel: Model comparison using the correlation matrix. The 
Bayes factor of the preferred model is shown in bold. 



tation through Euler angles a, 7 [? ] with < a < 27r, 
< P < TT and < 7 < 27r, realised by the Wigner 
D-functions [? ? ]. We have also let the amplitude of 
the power spectrum of initial perturbations vary slightly. 
Now the likelihood takes the form: 

1 1 T -1 

/:(A,17/e,a,A7) = /lo ^ I exp--[i?a] C^^ [Ra], 

where the A, ftk subscript denotes that the correlation 
matrix depends on the curvature and the amplitude of 
initial perturbations and R = R{a^ 7) is the rotation 
matrix (Wigner function). The amplitude A [? ] (the 
normalisation factor in CMBFAST with the option "un- 
norm") is related to the amplitude of fiuctuations at hori- 
zon crossing via |Ai^(/co)P = 2.95 x 10~^A. The data we 
use are the a^rnS extracted from the three-year WMAP 
ILC map whose noise at large scales is negligible. 

Marginalizing the likelihood function over the ampli- 
tude (with a fiat prior in the range [0.5, 1.2]) and the Eu- 
ler angles (we will explore the structure of the likelihood 
in more detail in [? ]), shows that the data favours spaces 
with little or no structure in their correlation matrices 
(Fig.[5j). Indeed, the likelihood function takes higher val- 
ues for less complex spaces, where correlations among 
various multipoles are weak. 

Finally, we calculate the Bayes factor for each topol- 
ogy using the same prior on Qk (Table [J). A compar- 
ison between the two panels shows that the inclusion 
of off-diagonal terms decreases the chances of all non- 
trivial topologies but the Quaternionic space. However, 
its Bayes factor is still too small to grant it any preference 
over the standard model. Changing to a Gaussian prior 
on Ho does not affect these results significantly. A com- 
putation of the evidence for the Truncated Cube space 
with n = 0.5 and n = 1.2 yields Bayes factors of 10~^ 
and 2.6 • 10~^ respectively. Given the effect of changing 
n, we expect this result to refiect mostly the change in 
the averaged diagonal power spectrum. 
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Figure 5: Likelihood based on the correlation matrix. The 
likelihood of the Poincare space is multiplied by 10^. 



CONCLUSION 

The non-trivial topologies we have examined are not a 
viable candidate for explaining the low-^ anomaly in the 
WMAP data. Bayesian model comparison shows that the 
flat, simply-connected ACDM model flts the data better. 
Considering only the power spectrum (assuming a con- 
stant n over our range of (3) topological spaces are penal- 
ized for requiring a low Hubble constant and tuning of 
the curvature, while a full analysis shows that the obser- 
vations do not support the extent of structure present in 
their correlation matrices. 

However, the possibility of a multi-connected topol- 
ogy is not necessarily excluded. We have limited our- 
selves to the simplest spherical multi-connected mani- 



folds; double- or linked- action manifolds have not been 
explored (however calculating their Laplacian eigen- 
modes could be computationally challenging). More im- 
portantly, our models required a low Hubble constant 
as a result of the geometrical degeneracy. It might be 
possible to adjust the power at large scales using other 
combinations of cosmological parameters, but this would 
require an additional mechanism to ensure that the shape 
of the power spectrum at small scales matches the obser- 
vations. Finally, it is possible that the primordial power 
spectrum is actually different from what we used. Lack- 
ing an established mechanism generating perturbations 
in closed spaces, we adopted the analogue of a scale- 
invariant power spectrum. If such a mechanism is ever 
conceived, it might predict a different form of a primor- 
dial power spectrum (although a change in the spectral 
index does not affect the pattern of correlations much). 
Another possibility would be to search for more informa- 
tion in other aspects of CMB observations, such as polar- 
isation, where a circles-in-the-sky approach might yield 
some results [? ]. But even in this case, a full analysis 
would have to follow the steps we described above. 

In conclusion, although the current data and the 
cosmological theories at hand do not support the case 
for a non-trivial topology, the number of possibilities left 
to explore suggests that the issue of the topology of the 
Universe is far from settled. 
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